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I.  Review  of  the  Project  Goals 

The  principal  purpose  of  the  project  as  specified  in  the  proposal  of 
September  1991  is  "to  evaluate  the  conceptual  design  of  an  acoustic 
system  which  is  purposefully  sensitive  to  the  identification  and  localization 
of  brief  transient  sounds". 

i 

In  the  Background  section  of  the  proposal  I  described  the  basic 
conceptual  elements  which  underlie  the  Silicon  Retina  (SR)  and  those 
which  underlie  passive  mammalian  hearing.  My  original  intention  was  to 
utilize  the  frequency  discrimination  capability  manifested  in  mammalian 
hearing  as  the  mechanism  on  which  to  base  an  acoustic  analogue  to  the 
SR  concept. 

Specifically,  mammals  have  sophisticated  frequency  dependent  sensors 
which  transmit  responses  to  discrete  and  well-resolved  frequencies 
upward  to  higher  processing  levels.  Each  intermediate  processor 
transmits  the  precise  frequency  information  upstream.  Thus  one  may  think 
of  acoustic  frequencies  as  being  analogous  to  visual  positional 
information,  since  in  the  ray  limit  of  optics  there  is  a  one  to  one 
correspondence  between  the  relative  locations  of  external  objects  and 
the  retinal  sensors  which  are  stimulated  at  any  given  time.  The  SR 
mechanism  was  to  be  implemented  by  comparing  individual  frequencies 
with  running  averages  over  groups  of  nearby  frequencies,  eliciting  a 
response  when  the  intensities  of  these  differed  by  a  prescribed  amount. 

The  next  step  in  the  program  was  to  remain  open-minded  regarding  the 
implementation  of  directional  acoustic  arrays,  with  the  idea  pf  steering 
an  array  making  use  of  the  "alert"  response  created  by  the  introduction 
of  a  non-background  signal.  In  this  way  it  was  envisioned  that  the 
"cocktail  party  effect"  would  be  achieved. 
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II.  Modifications  to  the  basic  strategy 

In  the  time  intervening  since  the  proposal  was  submitted  and  approved,  I 
have  become  aware  that  the  mammalian  auditory  system  is  far  from 
unique  in  purpose  and  function.  Evidently  hearing  has  been  "invented" 
somewhere  near  20  different  times  in  nature,  many  implementations 
being  radically,  rather  than  merely  superficially,  different  from  one 
another. 

Not  surprisingly  the  auditory  system  of  the  fishes  is  adapted  to  three 
dimensional  perception,  as  is  appropriate  for  the  underwater 
environment.  As  such  it  is  essentially  different  from  that  in  mammals.  There 
is  no  reason  of  principle  to  constrain  the  conceptual  framework  to  one  or 
another  biological  paradigm.  We  have  the  freedom  to  draw  desired 
elements  from  any  of  these.  In  what  follows  I  develop  reasons  why 
hearing  in  fishes  may  provide  methods  for  implementing  the  desired  “SR 
cocktail  party  effect  solver"  which  go  beyond  what  I  had  originally 
envisioned. 

III.  Physiological  and  physical  aspects  of  fish  auditory  systems 
A.  Physiological  characteristics 

Mammalian  hearing  achieves  its  precise  frequency  discrimination  through 
a  specialized  organ  known  as  the  Organ  of  Corti.  In  idealized  terms,  a 
shaped  nonlinear  membrane  has  evolved  which  maximally  responds  to 
different  frequencies  at  different  locations  along  its  length.  The  localized 
vibrations  of  this  membrane  stimulate  "hairs",  which  in  turn  stimulate  nerve 
endings  to  transmit  the  signal  up  to  higher  processing  levels. 

Fish,  it  seems,  have  no  direct  analogue  of  the  Organ  of  Corti,  nor  do  they 
possess  the  same  degree  of  frequency  sensitivity  and  discrimination 
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exhibited  by  mammals.  Their  needs  are  different.  To  begin  with  their 
body  mass  has  nearly  the  same  density  as  the  surrounding  acoustic 
medium,  and  in  consequence  they  may  be  roughly  described  as 
oscillating  bodily  in  response  to  the  acoustic  field.  Consequently  hearing 
is  possible  for  them  only  through  utilization  an  inertial  sensory  system. 

I  now  describe  this  inertial  system  in  broad  outline  [1  ].  Rather  than 
developing  a  nonlinear  frequency  discriminator,  the  fishes  evidently 
evolved  a  system  based  on  the  presence  of  a  "dead-weight"  bony 
structure  which  is  inertially  stable  in  the  presence  of  the  acoustic 
oscillation  of  the  bulk  of  the  fish.  As  in  mammals,  "hairs"  are  utilized  to 
sense  the  relative  motion  between  the  fish  and  the  bone.  These  hairs  are 
anchored  on  one  end  to  the  bone  and  on  the  other  to  a  membrane 
containing  nerve  endings.  Thus  one  may  think  of  the  tugging  of  the  hairs 
caused  by  the  relative  motion  of  the  endpoints  as  being  the  primary 
stimulus  to  “hearing". 

The  design  of  the  receptor  is  such  as  to  enable  directional  hearing.  Both 
the  magnitude  of  the  soundspeed  in  water  (many  times  that  of  sound  in 
air)  and  the  comparatively  small  distance  across  the  head  of  a  typical  fish 
render  the  mammalian  methods  of  directional  discrimination  (time  delay 
and  refraction)  ineffective.  Instead,  the  fish  (In  a  manner  of  speaking) 
hook  the  bone  to  hairs  pointing  in  different  directions.  This  wiring  scheme 
enables  an  x-y-z  decomposition  of  an  incoming  signal,  although  180° 
ambiguities  must  be  resolved  through  subsequent  processing.  (A  lucid 
discussion  of  a  model  for  how  th  is  is  accomplished  is  contained  in  Rogers 
et  al,  Ref.[l  ].  These  authors  also  discuss  frequency  discrimination.  I  will  not 
include  any  of  these  details  in  the  present  report.) 

A  moment’s  reflection  reveals  that  the  system  described  above  must  rely 
upon  pressure  gradients  rather  than  absolute  pressures  to  achieve  its 
function.  A  further  elaboration  of  the  scheme  utilizes  the  absolute 
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pressure  in  conjunction  with  the  three  Cartesian  components  of  the 
pressure  gradients.  While  the  details  of  the  biological  implementation  are 
interesting,  I  will  not  dwell  on  them  here.  Rather,  I  pass  to  a  discussion  of 
how  oceanic  noise  would  stimulate  a  system  designed  on  the  above 
conceptual  basis,  and  how  a  desired  signal  could  be  Identified. 

B.  Physical  considerations 

In  what  follows  I  will  assume  that  it  is  possible  to  construct  acoustic 
transducers  which  can  detect  the  absolute  pressure  level  and  the  three 
Cartesian  components  of  the  pressure  gradient  as  functions  of  time.  Such 
transducers  could  either  be  of  inertial  design,  mimicking  that  of  the  fish,  or 
could  utilize  (to  give  one  more  technologically  sophisticated  example) 
laser  interferometry.  Models  for  how  fish  might  process  acoustic  signals 
utilizing  such  "transducer"  systems  have  been  given  by  a  number  of 
authors.  Here  I  draw  upon  the  general  ideas  of  Rogers  et  al  [1  ],  although  I 
do  not  adhere  to  their  all  the  specifics  of  their  model,  since  my  aim  is 
different. 

For  the  purposes  of  the  present  discussion  I  also  limit  the  background 
noise  field  to  consist  of  surface  noise,  though  other  types  of  noise  will  be 
introduced  later  in  the  project .  I  do  not  foresee  that  their  introduction  will 
require  radical  modification  of  the  ideas  to  be  presented  here.  Again, 
many  researchers  have  contributed  to  our  knowledge  of  surface  noise, 
but  I  will  restrict  my  attention  to  the  formulation  of  Ffowcs  Williams  and 
Guo  [2.3]  in  the  present  work.  While  there  is  no  universal  agreement 
regarding  the  mechanisms  of  surface  noise  generation,  the  work  of  these 
authors  is  based  upon  the  Lighthill  definition  of  sound  sources: 

“the  source  of  sound  (is)  the  difference  between  the  exact  statements  of 
natural  laws  and  their  acoustical  approximations."^] 
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Based  upon  this  very  general  formulation  of  the  nature  of  the  sources  of 
transversely  homogeneous  and  isotropic  surface  noise,  it  can  be  shown 
that  the  far-field  pressure  in  the  water  p(z,x  ,t),  where  z  («.0)  is  the  depth 
and  x  is  a  two-dimensional  transverse  coordinate,  is  of  the  form 

i 

p(z,x  ,t)  =  x  (t-  |  z  |  /cw).  (l) 

Here  x  is  a  function  whose  argument  exhibits  retardation  of  the  arrival  at 
soundspeed  cw,  but  which  is  otherwise  arbitrary.  In  fact,  due  to  the 
irregular  motion  of  the  ocean  surface,  it  is  best  regarded  as  a  random 
function.  Proof  of  this  assertion  is  deferred  to  the  Appendix,  in  which  a 
more  general  situation  in  which  p  depends  on  x  is  also  discussed. 

It  follows  that 

vx  P  =  0; 

Vzp  =-cw-!  dtp.  (2) 

The  first  result  above  simply  reflects  that  the  surface  is  assumed  to  be 
transversely  homogeneous  and  isotropic,  at  least  in  a  statistical  sense,  so 
there  is  no  meaningful  reference  point  from  which  to  measure  x!  (If,  in  a 
different  context,  p  depends  on  x,  under  certain  conditions  the  second 
relation  above  can  be  generalized  to  Vx  p  s  ft  (-  cw_1  df  p)  in  the  far 
field,  and  the  ensuing  discussion  would  be  suitably  modified.)  Thus  at  this 
level  of  discussion  any  transverse  gradient  whatever  would  be  a  "signal” 
and  not  surface  noise.  In  practice  it  would  be  necessary  to  perform  a 
running  time  average  since  we  can  expect  x  dependent  fluctuations  to 
the  acoustic  pressure  to  occur. 
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For  convenience  let  us  now  define  <j>i  =  Vz  p  and  $2  =  -  cw_1  p. 

Construct  the  running  time  average 

1  T0  +  AT 

<4>1  >t  =  AT  Jdt  *1 

'  To 
1  T0  +  AT 

=  —  Jdt  at  <1)2 

To 

=  (<I>2(To  +  AT)  -  <j>2(To))-  (3) 

In  words  this  relation  states  that  the  running  time  average  over  the  z 
component  of  the  gradient  should  equal  a  specifically  weighted 
difference  of  the  absolute  pressure  measured  at  the  endpoints  of  the 
averaging  interval.  This  is  a  computation  eminently  suited  to  analog 
methods.  It  will  be  especially  useful  if  it  can  be  used  to  differentiate  signal 
from  noise. 

To  investigate  this  possibility  write  a  total  pressure  field  (suppressing  the  z 
dependence) 

P(t)  =  p(t)  +  S(t)  ©  (t  -t0)  0  (tf-t)  (4) 

where  p(t)  is  the  noise.  S(t)  is  a  signal,  and  ©  are  Heaviside  functions. 

Let  the  symbol  refer  to  the  total  field  P,  and  W  refer  to  S  using  the  same 
definitions  introduced  for  <|>  referring  to  p.  We  find 


7 


* 4,1  =  AT  ( +2(T°  +  AT)  *  4>2(to)  ) 

To  +  AT 

+  XT  Jdt  W]  (t)  ©  (t  -t0)  ©  (tf-t)  (5) 

To 

If  S  is  a  direct  signal  in  the  far-field  which  has  no  reverberation 
components  from  the  surface,  bottom,  or  intervening  structures,  it  will,  like 
satisfy  W]  2  In  this  case  the  integration  can  be  performed,  but  the 
endpoint  values  depend  upon  how  tQ  and  tf  are  related  to  the  endpoints 
of  the  averaging  interval.  Alternatively,  the  signal  may  not  satisfy  W]  2 
d\^2  Qt  Qll.  in  which  case  the  integral  does  not  simplify. 

For  definiteness  let  T0  refer  to  the  first  averaging  interval  in  which  the 
signal  appears,  T0  +  AT  >  tQ  >  T0  .  We  assume  further  that  at  earlier  times 
running  averages  were  performed  and  Equation  (3)  was  satisfied. 
Examining  Eq.(5)  we  recognize  that  4>2(To)  =  ^(To).  but  there  is  no  reason 
for  the  rest  of  the  right-hand  side  to  equal  $>2(Tq+  AT).  Thus  the  "test"  of 
Equation  (3)  based  upon  the  observed  total  field  fails. 

Let  us  next  suppose  that  the  signal  is  so  long  that  it  completely  overlaps  at 
least  one  subsequent  averaging  period.  In  this  case,  if  2  dt^2<  the 

right-hand  side  of  Eq.(5)  does  reconstruct  ~  (  <i>2(To  +  AT)  -  ^2(To) ) 

over  this  subsequent  interval,  and  the  "test”  is  a  success.  On  the  other 
hand  if  the  test  fails.  Thus  when  the  signal  totally  overlaps  the 

noise  there  is  no  chance  of  using  success  or  failure  of  the  test  in  order  to 
unambiguously  decide  whether  the  signal  is  present. 
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Finally,  it  should  by  now  be  clear  that  occurence  of  the  end  of  the  signal 
within  a  subsequent  averaging  interval  will  certainly  produce  a  failure, 
followed  by  a  success  due  to  the  ongoing  noise. 

To  summarize: 

•  Monitor  over  a  total  time  interval  much  longer  than  the  averaging 
interval  AT,  during  which  the  average  gradient  -to-  pressure  test  is 
succesful.  Provided  the  signal  does  not  accidently  begin  precisely  at 
the  start  of  the  averaging  step,  there  will  occur  a  first  averaging 
interval  in  which  the  test  fails,  indicating  the  presence  of  the  signal. 

•  If  the  signal  lasts  longer  than  successive  averaging  intervals,  the  test 

<  *1  >t  =  at  (  «2(To  +  AT)  •  ®2(To) )  wil1  succeed  over  those  intervals 

overlapped  completely  by  the  signal  provided  a  ^2  <  but  will  fall  if 
this  condition  is  not  satisfied. 

•  The  end  of  the  signal  sequence  will  be  marked  by  a  sure  fail  (barring 
happenstance)  followed  by  a  success  as  the  noise  persists. 

IV.  Remaining  tasks 

In  the  preceding  section  I  have  endeavored  to  report  on  new  directions 
in  which  the  sponsored  research  is  moving.  The  full  report  at  the 
termination  of  the  project  will  address  the  specific  Issues  raised  in  the 
project  proposal,  including  simulations  of  the  response  of  the  proposed 
transient  sound  selector  to  signals  superposed  on  realistic  models  of 
surface  noise.  In  performing  these  simulations  I  will  be  especially  focused 
on  issues  regarding  optimal  time-averaging  intervals  A,  given  the 
spectrum  of  the  noise,  and  on  sensitivity. 
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Appendix 

The  theory  of  Ffowcs  Williams  and  Guo  is  based  upon  partitioning  the 
active  region  of  sound  propagation  and  production  at  the  ocean  surface 
into  domains  z  >  A,  A  >  z  >  0,  and  z  <  0,  as  indicated  in  the  figure  below. 

In  the  first  and  last  domains  there  are  no  sound  sources,  so  the  free  wave 
equation  anpropriate  to  each  medium  obtains.  Let  p(z,  k;a>)  be  the 
Fourier  transform  of  the  pressure  field  with  respect  to  the  transverse 
directions  and  with  respect  to  time.  Thus  with  q  the  soundspeeds 


d2z  p  +  Yi2  P  =0; 


Yj2  =  (a)/Cj)2  -  k2. 

The  branches  of  the  Yi  are  chosen  so  that  the  pressure  vanishes  at  spatial 
infinity. 


Z  -  A 


l 


Acoustic  Reid 
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Since  we  are  dealing  with  fluids  in  motion  the  following  is  also  valid: 

Pw  df  V2  =  -  dz  p, 

where  V  is  the  fluid  velocity.  Thus  in  Fourier  space,  and  in  particular  at  z=0 


Yw  P(z=0)  =  -wpw  V2(z=0). 


Similarly 

Ya  P(z=A)  =  +  «  Pa  Vz(z=A). 


Using  these  relations  together  with  p  (z<0)  =  p  (0)  e-iYw  z  ,  and 
p  (z>A)  =  p  (A)ejYa  (z-A)  (  one  obtains  eventually 


p/zx  _  PaPwe'iYw  z 

;  PaYw  +  PwYa 


{^■p(0)-mvz(0)} 

Pa 


By  virtue  of  the  equation  before  this  one,  an  expression  more  symmetrical 
in  z=0  and  z  =  A  may  be  written,  but  this  is  unnecessary  for  what  follows. 

The  point  of  the  above  manipulations  has  been  to  express  the  acoustic 
field  within  the  water  in  terms  of  pressure  and  velocity  values  at  the 
surfaces  surrounding  the  acoustic  sources.  A  full  solution  to  the  problem 
would  require  examining  these  terms  in  detail.  But  the  general  character 
of  the  acoustic  field  is  already  embodied  in  the  result  above. 

To  see  this  in  brief,  we  follow  Guo  and  evaluate  the  inverse  Fourier 
transform  using  the  method  of  stationary  phase,  assuming  only  that  the 
expressions  involving  p(o)  or  p(A)  are  not  so  irregular  as  to  spoil  this 
estimate.  The  core  of  the  computation  is  In  the  evaluation  of  the  integral 
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l(x,y;z)  =  Jdk  exp  {i  k  •  (x  -  y)  -  z\Jv2  -  k2}  F(k(w) 


Here  F  is  a  regular  function  of  its  arguments,  and  the  structure  of  the 
phase  follows  from  the  definition  of  y  introduced  earlier,  v  =  aj/cw.  The 
phase  is  statignary  at 


k*  =  -v  (x  -  y)  /  R, 

where  R  =  yj  (x  -  y)2  +  z2.  Introducing  an  angle  tan  6  =  z  /  |  x  -  y  | ,  it  is 
possible  to  demonstrate  that  the  integral  takes  the  general  form 

I  a  f(0,(o)  exp(  -  i  wR/cw  )/R2 

where  f  is  polynomial  in  w.  Thus  the  inverse  Fourier  transform  back  into  the 
time  domain  will  introduce  Dirac  delta  functions  and  derivatives  of  these 
Delta  functions  with  argument  (t  +  R/cw  ),  i.e.,  the  "sources"  as  described 
by  p(0),  etc.,  are  retarded  by  the  time  needed  to  reach  x  from  y  with 
soundspeed  cw.  Schematically, 

p(x,z;t)  =  Jdy  9$>(*7y)]  0\  V<  t  +  R(x  -  y)/cw) 

The  form  of  <P  indicates  that  we  made  no  specific  assumptions  regarding 
the  transverse  (y)  dependence  of  the  sources.  There  is  no  general  reason 
for  the  spatial  gradients  to  be  simply  proportional  to  the  time  derivative.  It 
is  only  in  the  limit  of  homogeneous  and  isotropic  sources  that  we  can 
expect  the  y  integration  to  correspond  to  a  surface  averaging 
procedure,  leaving  the  result  utilized  in  the  main  part  of  this  report.  A 
similar  argument  can  be  made,  however,  for  directional  noise. 


